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Abstract. A general stochastic traffic cellular automaton (CA) model, which 
includes slow-to-start effect and driver's perspective, is proposed in this paper. It 
is shown that this model includes well known traffic CA models such as Nagel- 
Schreckenberg model, Quick-Start model, and Slow-to-Start model as specific cases. 
Fundamental diagrams of this new model clearly show metastable states around 
the critical density even when stochastic effect is present. We also obtain analytic 
expressions of the phase transition curve in phase diagrams by using approximate 
flow-density relations at boundaries. These phase transition curves are in excellent 
agreement with numerical results. 
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1. Introduction 

A traffic jam is one of the most serious issues in modern society. In Japan, the amount of 
financial loss due to traffic jam approximates 12 thousand billion yen per year according 
to Road Bureau, Ministry of Land, Infrastructure and Transport. 

Recently, investigations toward understanding traffic jam formation are done not 
only by engineers but also actively by physicists pQ. The typical examples are: 
developing realistic mathematical traffic models which give (with help of numerical 
simulations) various phenomena that reproducing empirical traffic flow [2], [31 HI [5J El [71 [8] ; 
fundamental studies of traffic models such as obtaining exact solutions, clarifying the 
structure of flow-density diagrams or phase diagrams [9j [10J [TlJ [12] ; analysis of traffic 
flow with bottleneck [131111] ; and studies of traffic flow in various road types [T5l [TEj [TT] . 

We can classify microscopic traffic model into two kinds; optimal velocity models 
and cellular automaton (CA) models. A merit of using a CA model is that, owing 
to its discreteness, it can be expressed in relatively simple rules even in the case of 
complex road geometry. Thus numerical simulations can be effectively performed and 
various structures of roads with multiple lanes can be easily incorporated into numerical 
simulations. However, the study of traffic CA model has relatively short history and 
we have not yet obtained the "best" traffic CA model which should be both realistic as 
well as simple. There are lots of traffic CA models proposed so far. For example 
Rule-184 [18] . which was originally presented by Wolfram as a part of Elementary 
CA, is the simplest traffic model. Fukui-Ishibashi (FI) model [|5] takes into account 
high speed effect of vehicles, Nagel-Schreckenberg (NS) model [2] deals with random 
braking effect, Quick-Start (QS) model [6J with driver's anticipation effect and Slow-to- 
Start (SIS) model [3] with inertia effect of cars. Asymmetric Simple Exclusion Process 
(ASEP) [ID] , which is a simple case of the NS model, has been often used to describe 
general nonequilibrium systems in low dimensions. An extension of the NS model, 
called VDR model, is considered by taking into account a kind of slow-to-start effect 
p. Recently Kerner et al proposed an elaborated CA model by taking into account 
the synchronized distance between cars [19]. Each of these models reproduces a part of 
features of empirical traffic flow. 

A fundamental diagram is usually used to examine whether a model is practical or 
not by comparing empirical data (see figure [1]) and simulation results. A fundamental 
diagram is called a flow-density graph in other words. And it chiefly consists of three 
parts: free-flow line, jamming line, and metastable branches. On the free-flow line, 
flow increases with density, while on the jamming line flow decreases with density. The 
critical density divides free-flow line from jamming line. Around the critical density, 
close to the maximal flow, there appear some metastable branches. In general, the 
metastable branches are unstable states where cars can run like free-flow state even 
if the density surpasses the critical density. Characteristics of individual fundamental 
diagrams for above CA models are as follows: the fundamental diagram of Rule-184 
has an isosceles triangular shape. The critical density in FI model, that divides the 
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whole region into free flow phase and jamming phase, is lower than that of Rule- 184, 
and the maximal flow reaches a higher value due to the high speed effect. The critical 
density and maximal flow in QS model are higher than those of Rule-184, because more 
than one cars can move simultaneously. The maximal flow of ASEP or NS model is 
lower than those of Rule-184 or FI model because of random braking effect. In the case 
of SIS model, a fundamental diagram is significantly different from these models. Its 
fundamental diagram shows metastable branches near the critical density [H H] . Because 
metastable branches are clearly observed in practical fundamental diagram (see figure 
[1]), we think realistic traffic models should reproduce such branches. We note that the 
study of these metastable branches in the high flow region is very important for a plan 
of ITS (Intelligent Transport Systems). The states in these branches are dangerous, 
because of the lack of headway between vehicles. On the other hand, they have the 
highest transportation efficiency because each vehicle's speed is very high regardless of 
the short headway. 

In this paper, we propose a general stochastic traffic CA model which includes, as 
special cases, well known models such as NS model, QS model, and SIS model. The 
paper is organized as follows. We first define our model in sec. 2. Next, fundamental 
diagrams, flow-a-/? diagrams and phase diagrams of this new model will be presented in 
sec. 3 and 4. Finally we show analytical expressions of phase transition curves in phase 
diagrams in sec. 5. 
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Figure 1. A fundamental diagram is plotted with use of empirical data. These 
data were collected at the Tomei expressway in Japan by Japan Highway Public 
Corporation. The road consists of two-lanes. Each points in the diagram corresponds 
to an average over a time interval of five minutes. The density is calculated from the 
flow divided by the average velocity. We can see metastable points which are inside 
the circle at the highest flux region of free flow. 
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2. A new stochastic CA model 

In 2004, a new deterministic traffic model which includes both of slow-to-start 
effects and driver's perspective (anticipation), was presented by Nishinari, Fukui, and 
Schadschneider [7]. We call this model Nishinari- Fukui-Schadschneider (NFS) model. 
In this paper we further extend NFS model by incorporating stochastic effects. 
First, the updating rules of NFS model are written as 
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in Lagrange representation. In these rules, x\ is a Lagrange variable that denotes the 
position of the ith car at time t. vf^ is velocity x\ — a;* -1 , and the parameter S represents 
the interaction horizon of drivers and is called an anticipation parameter. And we repeat 
and apply this rule again as < — v^f 1 in the next time. In this paper, we use a parallel 
update scheme where those rules are applied to all cars simultaneously. Rule (pQ) means 
acceleration and the maximum velocity is V max , (T5]) realizes slow-to-start effect, ([3]) 
means deceleration due to other cars, (jl]) guarantees avoidance of collision. Cars are 
moved according to ([5]). The characteristic of NFS model is the occurrence of metastable 
branch in a fundamental diagram because of the slow-to-start effect. It should be noted 
that the slow-to-start rule adopted in this paper ([2]) is different from the previously 
proposed ones [4J. 

Next, we explain a stochastic extension of NFS model. The model is described as 
follows: 
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We define independent three parameters, p, q, and r: the parameter p controls random- 
braking effect; that is, a vehicle decreases its velocity by 1 with the probability 1 — p. 
The parameter q denotes the probability that slow-to-start effect is on; the rule (jTJ) is 
effective with the probability q. r denotes the probability of anticipation; S = 2 with 
the probability r and 5=1 with 1 — r, hence the average of S is (S) = 1 + r. This 
value of S is used in rules (JTj) and (jHJ). We note that in order to reproduce empirical 
fundamental diagram the value of S should lie between 5 = 1 and S = 2 [7]. Thus we 
call this new extended model Stochastic(S)-NFS model. A remarkable point is that the 
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S-NFS model includes most of the other models when parameters are chosen specifically 
(see figure [2]). Here, "modified(m) FI model" means a modification of the original FI 
model so that the increase of car's velocity is at most unity per one time step, because, 
in the original model, cars can suddenly accelerate, for example, from velocity v = to 
v = Vmax, which seems to be unrealistic. 
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Figure 2. A reduction map of various traffic CA models. 



3. Fundamental diagram 

In this section, we consider fundamental diagrams of S-NFS model on periodic boundary 
condition (PBC). Figure [3] shows fundamental diagrams with maximum speed V max = 1. 
The parameter sets (a), (c), (g), and (i), correspond to Rule-184, SIS model, QS model 
with S = 2, and NFS model respectively. Figure H] shows those with maximam speed 
Vmax = 3. The case (a) is mFI model with Vmax = 3 and (i) is NFS model. 
The following two points are worthwhile to be mentioned; 

(1) Although free flow and jamming lines are both straight for r = or r = 1.0, the 
shape of jamming lines are roundish for r = 0.5 (see figures [3]- (d) or |H-(d),(e),(f)) 
like ASEP with random braking effect (p < 1.0). The origin of this behavior is as 
follows: Two successive cars can move simultaneously as long as S = 2. However, 
once S of a car moving rear changes to 1, this car must stop. Hence random change 
of S plays a similar role as random braking. 
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(d)q = 0.0, r = 0.5 



(e)q = 0.5, r = 0.5 



(f)9 
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Figure 3. Fundamental diagrams of S-NFS model with p = 1.0, L = 100 and 
Vmax = 1. We repeat simulations for 10 initial conditions at each density. The flow is 
averaged over the range from t = 50 to t = 100. When q ^ and for uniform initial 
distributions, metastable branches are clearly seen. 



(2) As is seen in figure Hr(h), metastable states "stably" exist even in the presence of 
stochastic effects, although in many other models metastable states soon vanish 
due to any perturbations. There are also some new branches around the critical 
density in our model. As seen in figured! these new branches may be related to the 
wide scattering area around the density 60-125 (vehicles/km) in the observed data. 
Thus we believe we have successfully reproduce the metastable branches observed 
"stably" in the empirical data. 



4. Phase diagram 

In this section, we consider phase diagrams of S-NFS model with open boundary 
condition (OBC) when V mauX = 1. Figure [5] shows our update rules for OBC: 

1) We put two cells at the sites —2, —1. At each site a car is injected with probability 
a. The car's velocity is set to 1. 
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Figure 4. Fundamental diagrams of S-NFS model with p = 1.0, L = 100 and 
V max = 3. Contrary to the cases with V mayi — 1, metastable states can remain even for 
random initial distributions. 



2) We put cells at the sites L, L + 1, and at each site a car is injected with probability 

1 — (3 (an outflow of a car from right end is obstructed by these cars). It's velocity 
is set to 0. 

3) Besides, we put cells at the sites L + 2, L + 3. At these two sites cars are always 

injected (this rule means that for cars at sites —2, . . . , L + 1, there exist at least 
two preceding cars at any time). It's velocity is set to 0. 

4) We apply updating rule of S-NFS to cars at sites —2, . . . , L+l. However slow-to-start 

effect ([TJ) apply only if < x\~ x and xl~] s < L - 1. 

5) We remove all cars at —2, —1 and L, . . . , L + 3 at the end of each time step. 

These rules are devised in order to avoid unnatural traffic jam caused by boundary 
conditions. Here, the rule 3) is needed in order that for cars at sites —2, . . . , L + 1 there 
always exist at least two preceding cars. We need to apply S-NFS update rules to cars 
at sites L, L + 1 because otherwise the velocity of the car at site L — 1 is determined 
irrepective of the state at L + 1 even when S = 2. 
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With use of this boundary condition, we can calculate flow-a-/3 diagrams (see figure 
[7]) form which phase diagrams are derived (see figures [T2|) . The phase diagram of ASEP 
has been already known (see figure E]). In this diagram, the whole a- (3 region is divided 
into two phases in the case p = 1.0: low-density (LD) phase which is "a controlled 
phase"; and high-density (HD) phase which is 11 f3 controlled phase". On the phase 
transition curve, the effect of a is balanced with that of /?, and then these two phases 
coexist. The flow-a-/3 diagrams of S-NFS model are considerably different form the 
result of ASEP (see figure Cj). The difference will be discussed in secEl We restrict 
ourselves to the case V max = 1 and p— 1.0 in the following. 



5. Phase transition curve 



In this section, let us derive an analytic expression of the first-order phase transition 
curve of S-NFS model by combining approximate flow-density relations at boundaries 
and fundamental diagrams. We extend the method proposed in articles [TTI [12]. The 
method consists of three steps: first, we relate the flow and density near the boundaries 
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Figure 5. Updating rules for an open boundary condition of S-NFS model when 
Vmax = 1. (a): Update rules (l6l)- (fTT|) are applied to cars within this region, (b): The 
rule of slow-to-start ([7]) is effective only when < x^ 1 and %l+g < L — 1. 
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Figure 6. Definition of flow-a-/3 diagram of ASEP (p = 1.0, q = r = 0.0) and a 
corresponding a-{3 phase diagram. The line of the first-order phase transition is given 
by a = 0. 
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(&)q = 0.0, r = 0.0 



(b)g = 0.5, r = 0.0 



[c)q = 1.0, r = 0.0 




(d)q = 0.0, r = 0.5 (e)q = 0.5, r = 0.5 (f)q = 1.0, r = 0.5 




(g)g = 0.0, r = 1.0 (h)q = 0.5, r = 1.0 (i)g = 1.0, r = 1.0 



Figure 7. Flow-a-/? diagrams of S-NFS model with p = 1.0 and Kn ax = 1- Note that 
flow suddenly rises near (3 = 1.0 in (b), (c), (h), and (i) because cars never take the 
slow-to-start mode even when q ^ 0. 
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of the system in OBC; next, in fundamental diagrams of PBC, we calculate the gradient 
of free-flow line and jamming line; finally we can obtain phase diagrams for OBC from 
the above results. 



5.1. Relation between flow and density near the boundaries 

First, we calculate each flow J\ and J r at left or right end of this system. Configurations 
in figure M contribute flow at each boundaries. Then with use of mean-field approach, 
the probability of the configration (a) is C(_i)[l — C( )]. Similarly the probabilities of (b), 
(c), and (d) are rc ( _i)C (0 )[l - c ( i)], c (L _i)[l -c (L) ], and rc^^c^l- c (L+1) } respectively. 
Here, denotes the probability of finding a car at the site j. Because C(_i) = a and 
C(l) = C(l+i) = 1 - (3, we get 

f J t = a(l - C( ))(l +rc( )) 
\ J r = /3[l + r(l-l3)]c {L - 1) 

where we assume C( ) = C(i). Note, we ignored possibility that a car at the site can 
not move by slow-to-start effect in (b). 
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Figure 8. Configurations of cars contributing to flow at left or right ends. Other 
configurations do not generate flow at both ends. 



5.2. Approximate expression of fundamental diagrams 

Next we approximate free-flow line and jamming line of PBC by straight lines as 

J f = p = C(p) (free-flow line) 

Jj = x(l — p) = x(l — C(l-i)) (jamming line) 

where Jf and Jj denote the flow of free-flow and jamming phases, and x is a magnitude 
of the gradient of jamming line. Here, we need to get the relations between x and q or 
r. The parameter x is related to the velocity of backward moving traffic clusters which 
we will calculate with use of mean-field approach. 

Figure shows the possible configurations at the right end of clusters. The symbols 
u g (t) g £ {0, 1, 2} mean the probability of each configurations at time t. The 
subscript g denotes the gradient of each configurations. For example, let us consider 
the probability that the configuration figure [9] (c) becomes (b) in the next time step, 
U2(t) — * u\{t + 1) (see figure [TOl) . In order that (b) occurs, the slow-to-start effect 
is not active for the car A (the probability 1 — q). The car B does not move when (i) 
the driver's perspective is not active (the probability 1 — r) or (ii) although the driver's 
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Figure 9. The connection between parameter u g {t) and configurations of cars. 
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Figure 10. The case that U2(t) becomes u±(t + 1). 



perspective is active, the car B still halts due to slow-to-start effect (the probability qr) . 
Thus we get the term U2(t)(l — q)(l — r + qr) which expresses the probability that u 2 (t) 
becomes Ui(t+ 1). Considering other possibilities in the same way, we get the following 
recursion relations 

u (t+l) = Ul (t)q(l-r)+u 2 (t)q (14) 
m(t + 1) = u (*)(l - r) + «i(t)(l - r + gr)(l - g + qr) + u 2 {t){l - g)(l - r + gr) (15) 
u 2 (t + 1) = u (t)r + ui(t)r(l - g)(l - g + gr) + u 2 (t)r(l - qf. (16) 



We solve the equations (fl4l . (IT3|) . and (TTB1) for the steady state tt g (t + 1) = u g (t) 
and we get 
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Thus, we take the expectation x = Ouq + lui + 2u 2 , and obtain the answer: 

1 + r — gr + g 2 r — 2g 2 r 2 
1 + q — qr + gr 2 — 2g 2 r 2 
The relations between x and g or r obtained from numerical simulations are plotted and 
compared with eq. ( 1201) in figure dU The theoretical curve ( 1201) gives good estimation of 
numerical results. 
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Figure 11. The relation between the gradient of jamming line x and probabilities q 
or r. The solid line denotes the approximate analytic result, and points denote results 
of numerical simulations. There seem some discrepancies in (b). The reason is that we 
assume the shape of jamming line in the fundamental diagram to be straight although 
it is actually roundish. 



5.3. Derivation of phase transition curve 

Finally, we derive phase transition curve for OBC from above results. From ( fl2l) and 
(j!3p . we obtain the following equations 

J f = Ji =^ C( ) = a[l - C( )][l +rc( )] (21) 
j J = j r z[l-c (L _i)]=/3[l + r(l-/3)]c(£_i), (22) 

which lead to 



or(r - 1) - 1 + Ja 2 (l + 2r + r 2 ) + 2a(l - r) + 1 
<*9 = ^ (23) 

C{L - l) = (3[l + r(l-(3)]+x (24) 

Using ( TTBl) and the equation J/ = Jj which is the characteristic of two phase coexistence, 
we can get C( ) = x[l — C(l-i)). Substituting into this equation, we can express (3 
by q, r, x, and a: 



_ 1 + r + y [c(o) - x] 2 (l + r) 2 + 4rx[c (0 ) - x]c (0) 

2r 2r[c(o) — x] 

where C( ) is written by a and r as equation (|23|) . This is the approximated curve of the 
phase boundary. 

It should be noted that (1231) and (1251) are a singular when r = 0.0. Expression (1251) 
becomes 

= ^TT^ ( 26 ) 

x(l + a) — a 

in this singular case. Figure [12] gives a comparison between the approximate phase 
transition curve ( 1251) and contour curves of flow- a-/? diagram obtained from numerical 
simulations. The phase transition curves are very close to those obtained from numerical 
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simulations. However we note that the shape of jamming line in the fundamental 
diagram is assumed to be straight as ([TBI . Therefore we see the approximate result 
is not good in the high flow region (where a and (3 are close to 1) for figure [12]- (d) 
because in this case the shape of jamming line is apparently roundish (see figure [3]- (d)). 

The phase diagram consists of LD and HD phases, but maximal-current (MC) 
phase does not appear. Driver's anticipation effect make the area of HD phase smaller 
whereas the slow-to-start effect gives the opposite tendency. Phase transition curve 
which is straight for ASEP is bent downward by driver's anticipation effect. However, 
if q 7^ 0, the slow-to-start effect becomes dominant as flow increases. 

The article [1_2] also discusses a phase diagram of a model which includes only 
slow-to-start effect. The structure of phase diagram is similar to those of figure [T2T(c). 
S-NFS model shows more variety of phase diagrams. 

According to p2] , for large (3 there appears the enhancement of the flow which looks 
like MC phase. The enhancement is due to the finite system size effect and vanishes 
with increasing system size. Figure [13] demonstrates the effects of system sizes on our 
results. The shape of flow does not change with increasing the system size. Since there 
does not appear a plateau at a = 0.75, we judge that MC phase does not exist when 
p = 1.0 for S-NFS model. 

6. Summary and conclusions 

In this paper, we have considered a stochastic extension of a traffic cellular automaton 
model recently proposed by one of the authors [7]. This model (which we call S-NFS 
model) contains three parameters which control random braking, slow-to-start, and 
driver's perspective effects. With special choice of these parameters, S-NFS model is 
reduced to previously known traffic CA models such as NS model, QS model, SIS model 
etc.. Hence S-NFS model can be considered as an unified model of these previously 
known models. Next, we investigated fundamental diagrams of S-NFS model. The shape 
of fundamental diagram of S-NFS model is similar to that with use of empirical data. 
Especially the metastable branches, which are indispensable to reproduce empirical 
fundamental diagrams, clearly seen even when stochastic effects are present. This 
robustness of metastable branches in this model is advantageous because empirical data 
shows metastability even though stochastic effects always exist in real traffic. Thus 
we can expect that S-NFS model captures essential feature of empirical traffic flow 
and is very useful for investigating various traffic phenomena such as jamming phase 
transition as well as for application in traffic engineering. Finally, we investigated phase 
diagrams of S-NFS model with Vm ax = 1 and p = 1.0. The analytic expression of phase 
transition curve is obtained from the approximate relations between flow and density 
near the boundary, together with approximate gradient of jamming lines in fundamental 
diagrams. The analytic phase transition curve successfully explains those of numerical 
simulations. Recently we have found that, in the case p < 1.0 MC phase appears in the 
S-NFS model. The analysis on the cases V max > 1 or p < 1.0 will be reported elsewhere 
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Figure 12. Approximate phase transition curves that are calculated by ([25]) . and 
contour curves of flow-a-/3 diagram (Fig. [7j). The theoretical curves are m good 
agreement with numerical results. 
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Figure 13. Flow is plotted against f3 at a = 0.25 and a — 0.75 when p = q = r = 1.0. 
The symbols + (•) represent data for L — 600 (L = 3000). Solid lines denote the 
approximate results (13]) with ([23]), (24]) and ([25]) . 
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